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We investigate improving the hadronic jet energy resolution using mass-constrained 

fits of tt" — >■ 77 decays using high granularity electromagnetic calorimeters. Single tt" 

studies have indicated a large potential for improvement in the energy resolution of 

TT^'s, typically reducing the average energy resolution by a factor of two for 4 GeV 

TT^'s. We apply this method to fully simulated multi-hadronic events with multiple 

TV 's with widely varying energies using the ILDOO detector model. Several methods 

flj . for identifying the correct pairings of photons with parent tt^'s were explored. The 

'■^ ' combinatorics become challenging as the number of tt^'s increases and we employ the 

' , Blossom V implementation of Edmonds' matching algorithm for handling this. For 

f— I ' events where both photons of the vr" are detected, the resulting solutions lead to an 

•1— ( ' improvement in the tt" component of the event energy resolution for 91.2 GeV Z" 

C/3 . events from 18.0%/y/E to 13.9%/VE using the ILDOO detector and its reconstruction 

O ' algorithms. This can be compared to a maximum potential improvement to 12.2%/y/E 

C/3 ' if all photon pairs are matched correctly using the current photon reconstruction. 



1 Introduction 



The concept of particle flow calorimetry [l][2] is widely seen as one of the most promising 
approaches for reconstructing the energy of hadronic jets in e~''e~ colliders. Much progress 
has been made on the particle- flow algorithms, particularly as described in [3]. The envisaged 
^■^^ ■ high granularity detectors open up new avenues to explore and further improve the event 

V/^ I reconstruction in the detector designs envisaged for future accelerator facilities such as the 

CN ' ILC and CLIC e+e" linear colliders. 

fT^ I The ability to reconstruct accurately individual photons in the midst of hadronic jets 

^^ i leads to the possibility to reconstruct individual di-photon resonances such as tt^'s. A typical 

Cn I multi-hadronic event consists on average of about 62% of the energy in charged hadrons, 

28% in photons, and 10% in neutral hadrons with large event-to-event fluctuations. The 
photons are primarily from ir^ — > 77 and we aim to reconstruct these intermediate tt^'s. 
The results from mass-constrained fits can significantly improve the prompt electromagnetic 
component of jet energy resolution. 



2 Overview 

The focus of this paper is to investigate the performance improvements that may be obtained 
for the reconstruction of the energy of multiple tt^'s in multi-hadronic decays of the Z'^ using 
mass-constrained fits based on a realistic full simulation and event reconstruction of the ILD 
detector. We first briefly discuss mass-constrained fits in section [3] and then review some 
of the results obtained with toy Monte Carlo studies for both single tt^'s and multiple 7r°'s 
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from Z" decays in sectional We then concentrate in section [5] on the resuhs obtained from 
friU simulation of the ILD detector for events with muhiple tt^'s from Z° decay at rest. A 
central issue is how to best pair up sibling photon candidates to parent tt'^'s. We use three 
pairing algorithms: 

1. Cheated (using MC truth information) 

2. A realistic matching algorithm with no recourse to truth information using Blossom V 

3. Same as 2, except that incorrect pairings in the solution are identified (using MC truth 
information) and replaced with the original photon measurements 

We examine the mass-constrained fitting performance for the three different pairing al- 
gorithms and compare the results with the standard calorimetric measurement with no 
mass-constrained fits. In order to make the study tractable with the current reconstruction 
software we have restricted to simulating only the prompt 7r° component of Z° decays (ig- 
noring potential complications from additional photon candidates originating from charged 
and neutral hadrons) and requiring that none of the photons convert in the tracking volume. 

3 Mass-Constrained Fits 

Given a 7r° — ^ 77 decay, we apply a mass-constrained fit where the invariant mass of the two 
photons is constrained to the tt" mass. The adjusted photon momenta result in an improved 
energy and direction resolution for the tt*^. In the fit we define the photon 3- vectors using 
the measured energy, E^ and spherical polar coordinates, 6 and (f) for the photon directions 
as estimated from the reconstructed cluster positions in the calorimeter and the assumption 
of the photons originating from the interaction point. The errors on the measurements are 
defined assuming that the errors on E^ 9 and are independent as is the case when the two 
photons are spatially resolved. The mass constraint relates the energies of the two photons. 
El and i?2, to the opening angle -012 between the two photons: 



m 



2 A T7. J7^ o;„2 



4Si^2sin^(^i2/2) (1) 



The fitting has been implemented with the MarlinKinfit package [4] . In addition to improv- 
ing the tt" energy resolution when the correct photons are paired the procedure leads to a 
well-defined error estimate for the fitted tt" energy which can potentially be used to better 
assign errors on the overall event reconstruction at the individual event level. However the 
latter has not been the focus of this paper. 

4 Toy Monte Carlo 

4.1 Single 7r° 

Single tt" studies have demonstrated strong potential improvements in t:^ energy resolu- 
tion [5] when the photon directions {9 and (p) can be inferred with high precision. To 
illustrate this we show below the dependence of the tt" energy resolution on the relevant fac- 
tors: photon angular resolution, tt" energy and cos^*, where 9* is the center-of-momentum 
(CM) decay angle. These particular estimates were obtained using a toy Monte Carlo which 
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assumed energy independent angular resolutions with the (j) uncertainty related to the 
uncertainty by Urj, = ag/ sm0. 

Table [T] shows as expected that the improvement in the stochastic parameter of the 
energy resolution, a where (Te/E = a/^/E, depends a great deal on the CM decay angle. 
The most dramatic improvement being for symmetric decay (cos 0* = 0) with almost no 
improvement at all for highly asymmetric decay. For each bin in cosO* the distributions 
arc rather Gaussian. For a reference case of a 4 GcV 7r° with a photon fractional energy 
resolution, cr£;/£' of 1Q%/\/E with E in GeV and a 9 resolution of 1 mrad, the 7r° energy 
resolutioro improves to 9.4%/v^ averaged over all CM decay angles while for symmetric 
decay the resolution improves to 4.2%/vi?. The precise measurement of the opening angle 
between the two photons, 012, is the reason for the improvement in energy resolution, and 
this can be influenced by improving the relative error on -012 - either by increasing the 
angle by decreasing the tt*^ energy for the same cos 9* or further improving the measurement 
precision of the photon directions. The behavior, averaged over CM decay angles, on tt'^ 
energy and angular resolution are also shown in the table. 



Dependence 
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Dependence 


on cos 9* 
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4 
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6.4 
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9.1 
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0.15 
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9.0 


1 


8.4 
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4.2 
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9.0 
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0.25 
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Table 1: Measured values of the stochastic energy resolution parameter, a, estimated from 
the rms of the distribution of the stochastically scaled energy residual, {E — Egcn) / \/Egen, 
where E is the fitted tt" energy estimate and Eg^n is the generator energy. The table 
shows the dependence on the angular resolution, tt" energy, and cos 9* respectively. Unless 
otherwise specified the table assumes a tt" energy of 4 GeV, ag = 1 mrad, and an intrinsic 
electromagnetic calorimeter resolution, adot set to 16%. 



4.2 Multiple n°'s: Z° -^ qq{q ^ u, d, s) 

We then explored the impact of tt*^ mass-constrained fits in Z° events with center of mass 
energy of 91.2 GeV with multiple tt^'s. The events were simulated in the same fashion 
as the single particle studies using a well behaved energy uncertainty following the usual 
asjE = a/y/E model and fixed angular resolutions. There are some important things to 
note about tt^'s in a Z" event with decay to light quarks. The average number of tt^'s is 
9.42 ± 0.32 [B] with a median energy of about 2.4 GeV. The energy distribution is highly 
skewed to lower energies with a few high energy particles to compensate. To see the impact 



''Estimated from the rms of tfic observed non-Gaussian distribution (consisting presumably of several 
Gaussians for different cos 9*) 
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on the electromagnetic contribution to the overall event uncertainty, we extracted the tt^ 
informatiorl^ from the events and considered the impact on the quantity a = aE/\E for 
7r°'s only. For these studies we only considered prompt tt^'s produced within 10 cm of the 
interaction point, and we also only considered tt^'s which decayed to two photons (98.8% 
branching fraction) . Several combinations of energy and angular uncertainties were studied 
and summarized in Table[2] The general trends are as expected. Improved angular resolution 
results in improved fitted energy resolution averaged over the randomized CM decay angles 
of the ensemble of 7r°'s. 



«dct (%) 


(j0 (mrad) 


afit (%) 


afit/^det 
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0.25 


4.94 


0.62 


16 


0.25 


9.75 


0.61 


32 


0.25 


21.2 


0.66 


8 


0.5 


5.31 


0.66 


16 


0.5 


10.2 


0.64 


32 


0.5 


21.7 


0.68 


8 


1.0 


5.78 


0.72 


16 


1.0 


10.9 


0.68 


32 


1.0 


22.5 


0.70 



Table 2: Energy resolution improvement for the prompt electromagnetic energy contribution 
to 91.2 GeV Z°'s. The resolution after fitting, agt, was estimated from the rms of the 
distribution of the stochastically scaled energy sum residual, {E — Egcn)/ \/Egcn- 



5 Full Simulation 

Full simulation testing of the fitting procedure was performed using the standard distribution 
of ilcsoft vOl-09 for the ILD detector version ILDOO described in the ILD Letter of Intent [7| 
with a 3.5 T solcnoidal field. This uses Mokka to simulate with GEANT4 the detailed shower 
development of photons in the material of the detector, and the Marlin based framework 
for reconstructing the corresponding energy deposits. The most relevant features of this 
detector are a hermetic electromagnetic calorimeter composed of Tungsten absorber plates 
and Silicon readout pads with cell sizes of 5 mm x 5 mm. The calorimeter has 29 longitudinal 
layers with fine sampling every 0.6 Xq for the front 20 layers and coarser 1.2 Xq sampling 
for the back 9 layers. The octagonal barrel part is located at a radius exceeding 1.85 m, and 
the endcaps are located 2.45 m longitudinally from the interaction point. 

Modifications to algorithms included revising the cluster position estimates to use the 
standard weighted center of gravity method. It is anticipated that further improvements 
to reconstruction of cluster position using shower fitting may lead to significantly improved 
angular resolution as discussed in [8|. 



'^For these studies we only retained prompt tt'^'s produced within 10 cm of the interaction point which 
decayed to two photons (rejecting the 1.2% of Dalitz decays). 
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5.1 Photon Resolution Modeling 

Performance of the fitting procedure relies on reliable estimation of the resolution and cor- 
rection of any biases. The fit variables are the energy and angular directions, therefore 
we modeled the resolutions for as, og, and a^. Estimates for resolution were obtained by 
fully simulating and reconstructing single photons using a broad range of energies (0.2, 0.4, 
0.8, 1.6, 3.2, 6.4, 12.8, 25.6 GeV). In order to be considered in estimating the resolution, 
we required that the single photons do not convert before the calorimeter and that there 
is only one reconstructed cluster. Energy resolution was set to depend on energy and be 
independent of angle: (Je/'/E = A + B\og2{E) with E in GcV. The end-cap and barrel 
regions were evaluated independently. For the barrel it was found that A = 0.1591 ± 0.0004 
and B = 0.0053 ± 0.0002 with a xVdof = 8.6/6. For the end-caps, A ^ 0.1468 ± 0.0004 and 
B = 0.0049 ± 0.0002 with a xVdof = 7.8/6. 

The angular resolutions, og and u,^, were determined separately in the barrel and end- 
cap regions taking into account their energy dependence. It was found that the angular 
resolution dependence on 9 could not quite be adequately described using simple low order 
polynomials. Therefore events were grouped into ten bins along the barrel, equally spaced 
in < \z\ < 2329 mm (we verified that the response is ^-symmetric). The end-caps were 
also grouped into ten bins equally spaced radially in the 500 mm < p < 1843 mm range. 
The energies used in this sample included the same energies as above, as well as 0.1 GeV. 
This sample space of E, p and E, z was then used to estimate the angular resolutions by 
performing bi-lincar interpolation between known values. 





Energy (GcV) 
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Inner End-Cap 


o-e 


0.1 


1.25 ±0.03 


0.87 ±0.03 


0.81 ±0.03 
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0.60 ±0.01 


(Tij, 


1.6 


0.80 ±0.01 


0.78 ±0.01 


0.80 ±0.02 


2.32 ±0.04 


ere 


25.6 


0.26 ±0.01 


0.32 ±0.01 


0.19 ±0.01 


0.31 ±0.03 


CT0 


25.6 


0.31 ±0.01 


0.29 ±0.01 


0.33 ±0.02 


1.04 ±0.03 



Table 3: Sub-sample of values for 9 and <j) resolutions for photons in ILDOO. Resolution 
units are mrad 



With these procedures the pull distributions of the measurements show that wc manage 
to model the energy and angular resolutions reasonably well. 

5.2 Multiple tt" Tests using tt^'s from Z*^ Decay 

The full simulation studies of Z'^ events arc based on simulating the multiple prompt-7r*^ 
component of 50,000 Z° events allowing the simulation to decay the tt^'s. We do not at 
present simulate the rest of the event (consisting mainly of charged hadrons and neutral 
hadrons). We require that photons are detected in a fiducial region of | cos 611 < 0.978 in 
order to be considered for fitting. In order to restrict to events where our current calorimeter 
based photon reconstruction estimates work, we also remove events where photon conversion 
(7 — >■ e+e~) or Dalitz decays (tt^ — ?► e^e^7) have occurred in the tracker by scanning for 
such decays using the Monte Carlo information. We expect that such events can be used in 
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the future. Additionally, there are regions of the detector where the reconstruction software 
does not record the Monte Carlo particle that deposited energy (e.g. regions near the beam 
pipe). Events where Monte Carlo particles could not be recovered are not considered in the 
analysis. 

5.3 Energy Sum Estimates 

In the full simulation of tt^'s there are complicating factors that must be dealt with. One 
of these is that especially low energy photons may not be detected or reconstructed by the 
software due to very low energy deposits or too few cells being hit. The overall reconstructed 
energy therefore usually underestimates the generator energy, and sometimes only one of 
the two photons from a tt" is detected. So we find it appropriate to use modified definitions 
of "generator energy" in some of our performance estimates. 

Symbol Definition < i^tot > < -E'tot' > 

Ego Energy from all tt'^'s in the event 21.2 18.1 

Egi Energy from reconstructed photons 20.7 17.7 

Eg2 Energy from tt^'s with both photons reconstructed 17.7 15.6 

EgO > Egl > Eg2 

Table 4: Definitions of generator energy sums used for analysis of Z^ events. Also given are 
the mean energy sums in GeV for Z^ events, < Etot >, and the same quantity after removal 
of conversions and Dalitz decays, < E^^^^' >. 



5.4 Using Truth Information 

First studies used Monte Carlo "truth" information to pair up photons properly. This 
should indicate the best possible performance of the mass-constrained fitting procedure. 
Identification of which Monte Carlo particle is associated with reconstructed particles from 
the software was done using energy deposits in the detector. Since overlap of energy deposits 
does occur, each cluster was associated with the Monte Carlo particle most responsible for 
the cluster according to energy. 

The process of matching daughter photons with parent tt^'s uses the Monte Carlo decay 
tree information. When a single Monte Carlo particle was assigned to more than one cluster 
(due to cluster splitting in the reconstruction software), the highest energy cluster was 
preferred for the matches. Additionally, to be consistent with how the matching problem is 
approached when not using truth information, after kinematic fits are performed, matches 
are cut if the fit probability is less than 1% and the energies of photons that arc not matched 
to their parent tt" is estimated using the standard calorimetric measurement. 

5.4.1 91.2 GeV Z" Performance 

The event selection yieldeqj 11,568 Z^ events at 91.2 GeV and standard reconstruction 
resulted in a Gaussian fitted arcco = 17.9 ± 0.1% for the prompt electromagnetic energy 



'^Notc that with an average prompt tc^ multiphcity of around 9.4, even with a relatively small probability 
per photon of converting in the detector material, the overall efficiency for no conversion is small 
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sum using generator energy Ego described above with a bias of (—7.3 ± 0.2%)/vi?- The 
ttreco value and the bias was estimated from a Gaussian fit within ±3 rms units of the 
observed distribution of the stochastically scaled energy residual, {E — Eg)/ y/E^. This 
same procedure is used for all subsequent resolution and bias estimates, with the fitted a 
giving the estimated resolution parameter, a, and the fitted mean giving the bias in units of 
I/VE. There is a significant bias in the measured energy due to missing photons using the 
EgQ estimator and so we have focusscd more on the Egi and Eg2 estimators in the following 
studies which give similar resolutions on measured energy but a much smaller bias of around 
{-1.0±0.2%)/VE. 

The matching process successfully reconstructed 75% of the generator tt*^ energy and 
improved the energy resolution to afit = 13.5 ± 0.2% (using Egi). This fraction is primarily 
due to missing low energy photons in the final reconstruction, which in turn results in lone 
unmatched photons. But the large impact of undetected photons becomes apparent if one 
uses Eg2 so that only tt'^'s with the potential for reconstruction are considered. This results 
in afit = 12.2 ± 0.2%. The improvement ratio afit/aroco = 0.68 is similar to the Z^ toy 
Monte-Carlo studies for angular resolutions of 0.5 and 1.0 mrad. 

5.5 Matching without Truth Information 

Reconstruction of tt'^'s without using truth information is done based on the x^ from the 
kinematic fits. The process involves five main steps: 

1. Assign all identified photons and the misidentified photons labelled "neutrons" as 
candidate photons. 

2. Perform mass-constrained fits on every pair of candidate photons and generate a list 
of candidate 7r°'s. 

3. Remove all candidate tt^'s where the fit probability is less than 1%. 

4. Select the "best" combination of tt'^'s such that each photon is used at most once. 

5. The energy estimate of candidate photons that are not assigned to any tt*^ in the 
solution reverts to the standard calorimetric estimate. 

5.5.1 Scoring with x^ 

There are many possible ways to score a combination of 7r° 's to determine which is "best" . 
The primary method we use is to sum the x^ from each fit in the solution to obtain an 
overall global fit probability for each solution. We then select the most probable solution 
that maximizes the total number of tt^'s reconstructed. Below is a table that illustrates this 
selection method: 

Soln Fits xVdof Prob 



a 


6 


5/6 = 0.83 


0.544 


b 


7 


8.2/7= 1.17 


0.315 


c 


7 


14/7 = 2 


0.051 



In this example, the maximum number of tt'^'s found was seven, so we eliminate option 
a. Between options b and c, the best x^ is b. 
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5.5.2 Combinatorics 

The number of potential solutions to consider after all candidate 7r° 's have been generated 
at step 2 above can be very large. For n photons with n being even, the number of solutions 
consisting of n/2 tt'^'s goes as (n — l)(n — 3)(n — 5)...l. Fortunately the 1% fit probability 
cut usually dramatically reduces this, but the issue of scale still exists. Even with just 10 
photons, typical of say a jet with 5 7r'''s, this number is 945. 

We have employed methods from graph theory to approach this problem. The problem 
of matching photons to parent pions can be modeled using a graph. Each vertex of the 
graph represents a photon and each edge represents the potential pairing between photons. 
We then assign a weight to each edge according to the chi-squared of the mass-constrained 
fit. This then becomes a problem of selecting edges such that we use each vertex (photon) 
at most once and we minimize the sum of the edge weights (min x^)- 

This problem is very close to an efficiently solved problem in graph theory classified as 
a weighted non-bipartite matching problem in which one tries to find a perfect matching 
which minimizes the sum of the edge weights. A perfect matching is one where each vertex 
is used exactly once. For a perfect matching to exist, the graph must have an even number 
of vertices. 

Clearly a real detector can not be expected to only detect an even number of photon 
candidates. Therefore, we use a method for converting any weighted graph into one where 
a perfect matching exists [5]. This involves creating a new graph by first duplicating the 
original, such that one now has two identical graphs. Then one creates new edges connecting 
each original vertex with its duplicate. This extra step allows photons to go unmatched 
which is a desirable feature given that we expect that some photons should not be used to 
reconstruct tt'^'s (e.g. due to missing sibling photons). 

Once the graph has been constructed, we can then use an efficient method for finding 
the solution. One such method, is the latest implementation of Edmonds' algorithm [10] 
called Blossom V, by Kolmogorov [Tl] which can solve the graph in polynomial time. The 
worst case performance for a graph with n vertices and m edges is 0{n^ m) but on average is 
much faster. This appears to be well suited to the problem as defined but is not necessarily 
convenient for systematically exploring alternative solutions. 

5.6 91.2 GeV Z" Performance 

As with using truth information, the standard reconstruction yielded 11,568 events and 
an energy resolution of arcco = 17.9 ± 0.1%. The above matching procedure was able to 
correctly reconstruct 51% of the generator tt^'s which represented 68% of the total 7r° energy 
per event {Egi). This resulted in an effective energy resolution of am = 14.6 ± 0.2% and a 
small bias of (-2.5 ± 0.2%)/V^. 

When one looks at 7r° energy where both photons are reconstructed by the software, the 
improvement to afit then becomes 13.9 ± 0.2% with a small bias of (—2.8 ± 0.2%)/\/E. In 
this case, 78% of the tt" energy that had the potential for fitting (i.e. Eg2) was fit correctly 
and afit/aroco = 0.77. 

5.7 Effects of Incorrect Photon Pairings 

The 7r° reconstruction process occasionally pairs photons incorrectly and thus generates 
TT^'s which have been fitted with an inappropriate mass constraint. In the above example, 
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am 
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EgO 


17.9 


-7.3 


15.4 


-8.9 


14.6 


-7.0 


13.9 -4.7 


Egl 


18.0 


-0.6 


14.6 


-2.5 


14.1 


-0.6 


13.5 1.8 


Eg2 


18.0 


-0.8 


13.9 


-2.8 


13.0 


-0.8 


12.2 1.8 



Table 5: Summary of mass-constrained fit procedure results in % applied to multiple tt^'s 
from Z'^ — ^ qq. Errors on the a and bias values are approximately 0.2% 



while 51% of tt^'s were reconstructed, an additional 19% were incorrectly identified. The 
effects of these errors can be seen by removing incorrect tt'^ 's and replacing with the original 
photons that were measured with the calorimeter and then comparing the final results. In 
the case of 91.2 GeV Z°'s these "corrected" events improve the effective energy resolution 
to 14.1 ±0.2% with a negligible bias (using Egi). When using the Eg2 estimator, and 
correcting for incorrect pairings, the resolution decreases to am = 13.0 ± 0.2% and the 
bias is also negligible. Overall this indicates that incorrect pairings do negatively influence 
the resolution performance in a significant manner and that there is some scope for further 
improvement if one can reduce the error rate. The bias from incorrect pairings is small and 
can likely be calibrated out. 



5.8 Results Summary 

A summary of the results comparing the resolution and bias estimates from the different 
estimators and pairing algorithms are given in Table [5] Individual distributions for the 
different pairing algorithms are shown in Figures 1-4 for the generator energy sum estimates 
which are corrected for undetected photons. It is clear from all the distributions that the 
a parameter determined from the Gaussian fits does not fully account for the observed 
peakedness of the distributions based on mass-constrained fits. Such behavior results from 
TT^'s with small | cos0*| being measured with much better energy resolution after the fit than 
the average fitted 7r°. 



5.9 Tuning the Algorithm 

There are three main parameters that currently can be tuned. Pmin is the minimum fit 
probability necessary to be considered part of a final match solution. The single photon 
chi-squared, Xy represents the edge weight between photons cloned during the graph modifi- 
cation phase that guarantees that a perfect matching exists. It essentially allows for photons 
to go unmatched in the final solution. E^™ is the minimum photon energy that is required 
for a photon to be considered available to be matched. In the case of 91.2 GeV Z°'s we 
found optimal values to be: Pmin = 0.01, Xj = 6.6348 (equivalent to 1% fit probability), 
and £'""" = 50 MeV. Any significant deviations from these values resulted in a worsening of 
performance. It is likely that one or more of these depends on characteristics of the events 
(e.g. number and energy of particles) and is an issue we are currently working on. 
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Figure 1: Measured stochastically scaled energy residuals [E—Eg)/ ^J'E^ for different pairing 
algorithms using Egi as the reference energy. The top panel shows standard reconstruction 
and the bottom panel shows the realistic matching algorithm. 
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Fitted Energy Residuals (Seen) (GeV) RIMS =0.1089 l\/ledian=0.0155 




(E^-EcO/VEgi 



Figure 2: Measured stochastically scaled energy residuals [E—Eg)/ ^J'E^ for different pairing 
algorithms using Egi as the reference energy. The top panel shows the cheated version and 
the bottom panel shows the realistic matching algorithm but with incorrect pairings replaced. 
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Measured Energy Residuals (Full PiO) (GeV) RMS =0.1445 Medlan=-0.0091 
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Figure 3: Measured stochastically scaled energy residuals [E—Eg)/ ^J'E^ for different pairing 
algorithms using Eg2 as the reference energy. The top panel shows standard reconstruction 
and the bottom panel shows the realistic matching algorithm. 
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Fitted Energy Residuals (Full PiO) (GeV) RMS =0.0990 Median=0.01 57 
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Figure 4: Measured stochastically scaled energy residuals [E—Eg)/ ^J'E^ for different pairing 
algorithms using Eg2 as the reference energy. The top panel shows the cheated version and 
the bottom panel shows the realistic matching algorithm but with incorrect pairings replaced. 
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6 Outlook 

Efficient reconstruction of fitted tt'^'s depends significantly on how well individual photons are 
reconstructed. There arc a number of areas where the photon reconstruction efficiency could 
be improved either by algorithmic improvements or potentially by improving the detector 
design. 

• At high energies, photon energy deposits tend to be split into several sub-clusters 

• At energies below about f 80 MeV the efficiency to reconstruct photons degrades reach- 
ing 50% at too MeV, and lower still for lower energies 

• Photons which convert in the tracker are currently not reconstructed well and should 
be amenable to reconstruction with high efficiency and excellent resolution 

The loss of low energy photons is of particular concern because it leads to cases where an 
energetic photon from a 7r° may more easily be paired with the wrong photon. Cases where 
one of the photons from a 7r° converts in the tracker and Dalitz decays should lead to a 
significant improvement in the overall energy resolution on the tt*' through mass-constrained 
fitting. Other directions for improvement include improving upon the photon position re- 
construction, the use of the information from the fitted uncertainties in assigning jet-specific 
energy resolution uncertainties, and the assessment of performance for events with higher 
energy and particle multiplicity. 

7 Conclusions 

We have shown that with a full simulation of the ILD detector response to multiple tt^'s 
that it is feasible to improve substantially the prompt electromagnetic component of jet 
energy resolution using mass-constrained fits. For events where both photons of the tt" 
are detected, the resulting solutions lead to an improvement in the tt*^ component of the 
event energy resolution for 91.2 GeV Z^ events from 18.0%/^/E to 13.9%/a/E using the 
ILDOO detector and its reconstruction algorithms. This can be compared to a maximum 
potential improvement to 12.2%/\/E if all photon pairs are matched correctly using the 
current photon reconstruction. 
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